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Abstract 

We investigate the presence of quantum chaos in the spectrum of the bidimen- 
sional Fermi hquid by means of analytical and numerical methods. This model 
is integrable in a certain limit by bosonization of the Fermi surface. We study 
the effect on the level statisticsof the momentum cutoff A present in the bidimen- 
sional bosonization procedure. We first analyse the level spacing statistics in the 
A-restricted Hilbert space in one dimension. With g2 and (74 interactions, the level 
statistics are found to be Poissonian at low energies, and G.O.E. at higher energies, 
for a given cut-off A. In order to study this cross-over, a finite temperature is intro- 
duced as a way of focussing, for a large inverse temperature on the low energy 
many-body states, and driving the statistics from G.O.E. to Poissonian. As far as 
two dimensions are concerned, we diagonalize the Fermi liquid Hamiltonian with 
a small number of orbitals. The level spacing statistics are found to be Poissonian 
in the A-restricted Hilbert space, provided the diagonal elements are of the same 
order of magnitude as the off-diagonal matrix elements of the Hamiltonian. 
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1 Introduction 

The ideas of quantum chaos have recently been applied to the field of strongly correlated 
electron theory [jl| . These methods allow a non perturbative description of the statistical 
properties of a Hamiltonian of strongly correlated electrons, and may be a tool to extract 
some information from finite size systems. The aim of this article is to understand whether 
the methods of quantum chaos can shed some new light on the problems of strongly correlated 
bidimensional Fermi systems. The question whether the bidimensional Hubbard model is a 
Fermi liquid or not is controversial. From a theoretical point of view, Anderson [^] suggests 
that the ground state of the bidimensional Hubbard model is similar to the one of the one- 
dimensional Luttinger liquid, with Z = in the thermodynamic limit, whereas Engelbrecht 
and Randeria argue that the Fermi-liquid theory is not violated, so that this question is 
controversial Numerical computations lead also to controversial answers. For instance, 
Dagotto et al [0] |^ observe quasiparticules in the bidimensional t-J model, whereas Sorella 
emphasizes the Luttinger liquid behaviour. From the point of view of quantum chaos, one 
has to answer the question: do level spacing statistics contain information about the Fermi 



1 



liquid behaviour or non Fermi liquid behaviour of the t-J model? Before analysing the level 
spacing statistics of the t-J model we analyze in this paper models with a well established 



physical content, that is the Fermi liquid. As we shall see, we can answer only partially the 
previous question: one can detect integrable modes at low energy, but quantum chaos by itself 
does not give information as to whether these modes are quasiparticule modes or Luttinger 
liquid modes. However, there is one case in which one can conclude from quantum chaos: in 
the absence of integrable degrees of freedom at low energy, one can conclude to the absence of 
a Fermi liquid at low energy. Notice here the difference between the approach of [|l| and our 
point of view. In all the energy levels of the t-J model are analyzed on an equal footing, 
whereas we shall focus essentially on the low energy degrees of freedom of the Fermi liquid. 



We have established in the first article of this series |jTT|], that the level spacing statistics of 
the bidimensional Landau Hamiltonian are poissonian. In this model, the quasiparticules are 
in a one to one correspondence with the non interacting gas of spinless electron excitations, 
and thus occupy orbitals labelled by the same quantum numbers k as in the case of the gas. 
The quasiparticules interact among themselves in a diagonal manner 

H[{Sn]^}] =J2^^^'^^ + H /kk'^nk^^k'. (1) 

k (k,k'> 

A numerical computation proved that the level spacing statistics of the Hamiltonian (0) are 
poissonian in two dimensions, and close to Poisson in one dimension, with our truncation of the 
Hilbert space. The notion of level spacing statistics seems to be relevant for quantum fluids, 
and one has to distinguish between one and two dimensions. The link between the breakdown 
of the Fermi liquid picture in one dimension and the level spacing statistics has already been 
studied in reference |]13[. As far as two dimensions are concerned, we showed in the paper I 



111 , that the Landau liquid was characterized by its generical integrability, namely by Poisson 
level spacing statistics in two dimensions. Again, the level spacing statistics are a good tool 
to see whether the Fermi system is at the Landau fixed point or not. For the liquid to be 
a Fermi liquid, one should be able to generate quasiparticules of the interacting system by a 
switching on procedure. In the framework of the Landau theory, the success of the switching- 



on procedure |T^ suggests the conservation of the number of conserved quantities at low 
energy during the switching-on procedure, namely that the level spacing statistics of the gas 
and of the interacting liquid belong to the same universality class. Typically, adding a non 
diagonal perturbation to (|l]) and obtaining a Gaussian Orthogonal Ensemble (G.O.E.) level 
spacing statistics would mean a departure from the Landau fixed point. To be at the Landau 
fixed point, one must exclude strong correlations between levels. The corresponding generic 
level spacing statistics are poissonian. However, the case of the one dimensional Landau 
Hamiltonian shows that the statistics may not be exactly poissonian, but close to a Poisson 
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law. Poisson or close to Poisson level spacing statistics are not a sufficient condition for 
the system to exhibit quasiparticules since one could imagine a situation in which the same 
mechanism as in one dimension for the breakdown of the Fermi liquid holds, so that the level 
spacing statistics would remain poissonian even in the non Fermi liquid case. However, the 



adiabatic procedure, as emphasized by Anderson |Tj], is performed within a finite characteristic 
time 1/e, and the thermodynamic limit is taken for a finite e which is then much larger than 
the typical level spacing. So the adiabatic continuation does not generate true eigenstates (for 
this to be the case, one should take the limit e ^ ffist, before the thermodynamic limit). 
This is why the notion of adiabatic continuation advocated to introduce Fermi liquid is weaker 
than a similar requirement for all the eigenstates taken separately. So, this is why we should 
ask whether a Fermi liquid follows Poisson or G.O.E. statistics. 

The aim of the present paper is to study the level spacing statistics of the Hamiltonian of 
spinless electrons 

H = H^ + H^ (2) 
H' = E<k)c+Ck (3) 

k 

= ^ E /k,k',q4+q4'-qCk'Ck, (4) 
k,k',q 

where V is the volume of the system. The Hamiltonian (|I|) describes a Fermi liquid for 
times inferior as the decay time of the quasiparticules, and leads to good t her mo dynamical 
predictions However, the phenomenon of decay of the quasiparticules is not described by 



the Landau form ([l|). Since quasiparticules are true eigenmodes of (P, they have an infinite 
life-time. By contrast, the form takes into account the decay of the quasiparticules and it 
is indeed possible to calculate the decay rate of the quasiparticules from the Hamiltonian (Q) 



r^|. Even though the Hamiltonian it is not diagonal, it can be brought to diagonal form 



using some asumptions and by bosonizing the Fermi surface |T8[. We shall review the main 
assumptions and the bosonization of (Q) in the thermodynamic limit and in the limit of a zero 
curvature of the Fermi surface. One of the ingredients of the solution via bosonization is the 
existence of a momentum A which determines the single particule states which participate 
in the formation of bosons in the vicinity of the Fermi surface. A is a necessary ingredient, 
because of the curvature of the Fermi surface, but is also a source of difficulty since the objects 
with a true bosonic character in the limit A +00 are no more exactly bosonic in the limit 
of a finite momentum cutoff A. 

Of course, numerical computations cannot be performed in the thermodynamic limit since 
one can only diagonalize matrices of size 2000 by a Jacobi method. This technical limitation 
imposes to work with small systems and to impose a drastic cutoff A, so that the Hamiltonian 
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(|^) is no more integrable in the framework of the Hilbert space of the numerical diagonal- 
isations. In the case of the presence of a cutoff, and for finite size systems, one may thus 
expect the universal Poisson level spacing statistics to be replaced by the universal G.O.E. 
level spacing statistics of the form 

P{s) = -sexpi^-—j. (5) 

It is necessary to determine the importance of the cutoff A on the level spacing statistics and 
whether the statistics evolve towards G.O.E. level spacing statistics as A is reduced. This has 
not only a numerical interest, but also a physical one. For instance, lattice models have a 
natural cutoff A ~ 1/a, where a is the lattice spacing. To answer this question, we come back 
to one dimensional spinless fermionic systems, but with a cutoff A and study the level spacing 
statistics for g2 and (74 interactions. As we shall see, the level spacing statistics are drastically 
affected by the presence of the cutoff since the bosonic superpositions of particule-hole pairs 
which guarantee the integrability, do not survive in the presence of the cutoff. If the cutoff A 
is fixed and the size L of the system increases, one expects for a crossover from a G.O.E. level 
spacing statistics in the high energy part of the spectrum to a Poisson level spacing statistics 
at low energy, which shows that finite size effects are drastic in the presence of the cutoff. In 
order to determine the crossover scale between the Poisson regime and the G.O.E. regime, we 
introduce temperature-dependent level spacing statistics. The statistics at finite temperatures 
are found to pass from a non integrable statistics at high temperatures to Poisson statistics 
at low temperature. 

As far as two dimensionnal systems are concerned, we carry out numerical diagonalisations 
of a small system of electrons, with a small number of orbitals. In these conditions, the 
Hamiltonian is non integrable. We found Poisson level spacing statistics even for a small 
system. We attribute this property to the fact that the interactions generate extra diagonal 
matrix elements in the Hamiltonian which are in competition with diagonal interaction matrix 
elements, but the interaction connects only a small number of states, leaving a lot of zero 
matrix elements in the Hamiltonian. Provided the off-diagonal matrix elements are of the same 
order than the diagonal ones, the statistics are dominated by the diagonal matrix elements. 
However, by keeping only off-diagonal interaction matrix elements, we were able to exhibit 
G.O.E. level spacing statistics. 

This paper is organized as follows. We first treat the one dimensional case, in the Hilbert 
space restricted by the momentum cut-off A. The results established in one dimension are 
helpful in two dimensions since the presence of the momentum cut-off is also a source of non 
integrability. In a second step, we come to the bidimensional Fermi liquids, with the study of 
the level spacing statistics of the Hamiltonian (|^). 
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2 Level spacing statistics of 1 D spinless fermion sys- 
tems in the A-restricted Hilbert space 

2.1 Bosonization and level spacing statistics in the unrestricted 
Hilbert space (Luttinger liquid) 



Using the bosonization procedure of Haldane |2^, we can solve the one dimensional spinless 
Luttinger liquid with a linear dispersion relation and g2 and interactions. The Hamiltonian 
reads, in term of fermions 

H = Vf "^{ak - kp) : claCk,a ■ lid-igSaa' + g2qSa-a') Pq,aP-q,a' , (6) 

k,a a,a' 

where the label a indexes the branch a=R(ight),L(eft)= +1, —1. The interaction g^ describes 
the diffusion of two fermions on the same right or left branch whereas g2 describes the diffusion 
of two fermions belonging to the right and left branch. The density operators are defined as 

Pq,a = ^ • '^k+q,a'^k,a (7) 
k 

and obey bosonic commutation relations 

[Pq,a, Pg'^a'] = -Sa,a' 5q,q' , (8) 



which allows the definition of boson operators 

a; = (^)-"^^(ag)pg,„, (9) 



(^^1/2- 

L\q\ 



and make it possible to diagonalize the Hamiltonian (|^) via a Bogoliubov transformation: 

H = Eo + Y^iK^ + T^K^' + ^JJ')^ (10) 

q 

where 

6^ = cosh (yjgtt^ — Sinh (y9qa_g (11) 

tanh2(/?„ = ^ — (12) 

Vf + giq 

uJq = \{vf + g,,y - {g2q)V\<l\ (13) 

vs = {{vf + 9ao)' ~ {92o?y^' (14) 

vn = vsexp-2Lpo (15) 

vj = vsexp2Lpo. (16) 
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The interaction functions (72(0') and g^iq) are supposed to tend to a constant in the hmit g — 0, 
and to zero in the hmit q +00. Their decrease is controlled by the impulsion scale 271 /R, 
where i? is a given length scale. In the thermodynamic limit {L ^ R), one can define an 
effective low energy Hilbert space. To do so, we assume that Uq ~ ujq for all the wave vectors 
q. Of course, this approximation is only valid if the temperature is low enough. In this limit, 
the Hamiltonian (0) depends only on two velocities vn and vj: 

H = Eo + Y.ivNVjf'Wh'^qh, + ^{vnN^ + vjJ^). (17) 



9^0 



In one dimension, the link between the breakdown of the Fermi liquid picture and the level 
spacing spacing statistics is well understood. The absence of a Fermi surface in one dimension, 
in the presence of long range g2 and 5^4 interactions was known to Dzyaloshinskii and Larkin 
in the 70 's |T^. The breakdown of the Fermi liquid picture is already present at the level of 



the static correlation functions, and in the Hilbert space of low energy. The usual infrared 
divergencies are governed by the g — limit of g2{q) and ^'4(5), and are related to the 
orthogonality catastrophy and the absence of a Fermi surface. The infrared spectrum is 
described by the Hamiltonian (|T^). If the total charge and current quantum numbers are given, 
the q ^ ^ excitations are bosons with a linear dispersion relation, but with a renormalized 
velocity. The departure from the gas behaviour is measured by the anomalous exponents 
which appear in the static Green's functions. The interaction energy scale associated to this 
static breakdown of the Fermi liquid is given by [p 



However, as emphasized in [T^, the static part of the Green's function does not contain all 



the physics of the breakdown of the Fermi liquid theory. From a calculation of the two points 
Green's function, and from a calculation of the switching-on pocedure, we found in [13| a 
dynamical breakdown of the Fermi liquid picture, controlled by the interaction scale 

qiv- ^ (19) 

where a is defined by 

g2g = g2oil-iqRr)- (20) 

The static process of breakdown of the Fermi liquid picture does not modifie the level spacing 
statistics, which remain singular, as in the case of the gas, since the dispersion relation of the 
bosons remain linear. By contrast, the dynamical breakdown of the Fermi liquid appears to be 
related to a cross-over in the energy level spacing statistics, namely from singular level spacing 
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statistics to generical Poisson level spacing statistics. As far as symetries are concerned, the 
gas case and the Luttinger liquid possess the conformal symetry. The static breakdown of the 
Fermi liquid preserves the conformal invariance. The theory is the antiperiodic-antiperiodic 
sector of the compactified boson, with an interaction-dependent radius of compactification. 
The dynamical breakdown of the Luttinger liquid corresponds to a massless loss of conformal 



mvariance 



2.2 Fourier transform of the gas spectrum and temperature-dependent 
level spacing statistics 

We first wish to characterize the spectrum of the gas in the absence of interactions, a momen- 
tum cutoff A and a finite size L. The idea is to find a criterium to detect when a boson of given 
wave vector q is present or not in the system. Since g2 and (74 interactions are diagonal on the 
basis of bosonic excitations, this tool is a good way to characterize the degree of integrability 
of the model with a momentum cut-off. We consider a one branch model, and the number of 
quantum states is simply the integer part of 2AL/27r. The spectrum is made up of levels with 
an equidistant separation 2'KVp/L. In order to characterize the spectrum, we use its Fourier 
transform 

hir) = / Y: ^(^ -E.) = Y: e (21) 



In the thermodynamic limit. 



{E^} {Ei} 



hoo 



/oo(r)=/ e^--E'^(^-^/ E (22) 

J —00 r 1 „ _n 



leading to 



This function presents poles which are characteristic of the bosonic modes for r„ g = 2iTn/vfq, 
with n an integer. If we rescale r by a factor 27r and choose Vf = 1, we obtain a pole for 
Tn,g = n/q, that is for every rational number. We plotted |/oo(t)| for 5 bosonic modes on figure 
|l[ One can easily recognize the different bosons in the sequence of poles. In the presence of 
interactions, the poles are deplaced. We also studied the Fourier transform of the spectrum 
in the case of a finite size system. The result is depicted on figure ^ We can see that the 
poles are not so well defined. However, we can recognize the formation of peaks which replace 
the poles, and attribute well defined boson wave vectors to some peaks. The corresponding 
truncated density operators are defined by 

PA(g) = E^a(^)^a(A; + q)4+gCk, (24) 
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where 9\{k) = 6 {A — \k — kf\). We can ask under which condition a peak corresponding to 
a truncated boson Pa(q') appears on the modulus of the Fourier transform of the spectrum 
\f 1(^)1 ■ To answer this question, we assume that the number Ug of bosonic excitations of wave 
vector q appearing in the presence of a cutoff is such as qug ~ 2A, so that fL{ji,q) ~ 2A/g for 
the peak at r = ri^g. A boson of wave vector q is well defined provided fL{Ti,q) ^ 1, that is 
if A ^ q/2. For a given cut-off, the statistics are expected to be Poisson at low energies, and 
G.O.E. at higher energies. This is due to the fact that truncated bosons are created essentially 
at low energy. In order to test this idea, we introduce temperature-dependent level spacing 
statistics, which are defined as follows. If {Ei} is the full spectrum and {cj} are the energy 
levels after the smoothing procedure |2^, the density of level spacing P{s) is 




^/3(^)=(Eexp(- r^V^1 l Eexp(- /^^\+^M ^(.-(6.+i-6.)). (25) 



Moreover, for the statistics to be comparable, one needs to scale P/3{s) and s such as {Pp^s)) = 
1, and (sP/3(s)) = 1. The reason why we have to rescale these quantities is that they were 
equal to unity after the smoothing procedure, which was carried out in the zero temperature 



limit. This property is no more valid with the statistical weights of equation (^51) . For an 
infinite temperature {jS = 0), we recover the usual level spacing statistics. As the inverse 
temperature f3 increases, the low energy levels carry more and more statistical weight. As we 
saw from the description in terms of the truncated bosons (0), the spectrum is expected to 
be integrable at low energies, so that the level spacing statistics should evolve from G.O.E. 
statistics, or at least intermediate statistics (that is, with < P(0) < 1), to Poisson level 
spacing statistics as f3 increases. 

2.3 Level spacing statistics with §2 interactions in the A-restricted 
Hilbert space 

The numerical method to compute the level spacing statistics of two branch models with a 
momentum cutoff A and g2 interactions consists in tabulating the states of the Hilbert space 
by generating all the different fillings of the two branches, the number of particules on each 
branch being kept constant. The second step consists in computing all the matrix elements 
of the hamiltonian and diagonalising the Hamiltonian by the use of the Jacobi method. The 
interactions are of the form 

qL 



92^ = 92oew-[^). (26) 

The scale of the interactions R is chosen equal to A, so that all the matrix elements of the 
interaction are important. The evolution of the spectrum as a function of g2o is plotted on 



figure a. This figure is to be compared with the figure 2 of reference [13|, where we have plotted 
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the evolution of the energy levels as a function of g2o but in the absence of cutoff. The obvious 
difference between the two plots is the presence of level repulsion with a cutoff A, and the 
existence of level crossings in the absence of cutoff. The level spacing statistics corresponding 
to a spectrum analog to the one of figure |^ are plotted on figure |[ They are in good agreement 
with G.O.E. level spacing statistics. In our computation, KL/2'k = 3, so that the criterium 
A ^ 27r/L is not verified. We could not go to higher values of AL/2'k because increasing the 
ratio AL/2iT increases the size of the Hilbert space and numerical diagonalizations are no more 
possible. We conclude that the statistics depend drastically on the length of the system, if the 
momentum cutoff is fixed. If we apply the ideas of temperature-dependent statistics of level 
spacings to the system, we find that the statistics of level spacings are driven from a G.O.E. 
law at /? = to a Poisson law as (3 increases, as depicted on figure ^. This evolution reveals 
the fact that "nearly integrable" degrees of freedom exist at low energy. These bosonic degrees 
of freedom are exactly integrable in the absence of cut-off. In term of classical trajectories 
(provided one is able to find a classical phase space for the Fermi hquid !), this situation 
corresponds to the existence of conserved torii at low energy, which transform into chaotic 
trajectories as the energy increases. The cross-over temperature scale shall be derived later. 

2.4 Level spacing statistics with §4 interactions in the A-restricted 
Hilbert space 

We now consider the case of (74 interactions. The interaction Hamiltonian is 

=Y1 94qCt+q4'-gCk'Ck. (27) 

The interest of the g4 term is that we can use only a one branch model, and we can reach 
higher values of the ratio AL/2tt without increasing the size of the Hilbert space. We could 
reach AL/2tt = 12 in a sector of total momentum P = 24.27r/L. P is the impulsion with 
respect to the fundamental. As long as P < AL/vr, one has generated the complete Hilbert 
space in the absence of interactions. This fact motivates the choice P = AL/tt. The condition 
AL/tt = 24 ^ 1 is respected. However, we did not find Poisson statistics, but statistics 
which are intermediate between a Poisson law and G.O.E. statistics. Namely, the value of 
the density of normalized zero crossings is not 1 but 0.5. The statistics are plotted on figure 
^ This computation gives an idea of the extension of the crossover as a function of AL/2'7r, 
since we obtain intermediate statistics for AL/2tt = 24. We computed also the temperature- 
dependent level spacing statistics. We found a cross-over to Poisson statistics as the inverse 
temperature increases, as plotted on figure 0. This cross-over is much more rapid than in the 
case of the g2 interaction. However, the values of the cutoff in the g2 case and the g^ case are 
not comparable. 
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2.5 Cross-over scales 



In order to get an idea of the temperature cross-over between the G.O.E. and Poisson regime, 
we look for the energy scale fc^T* below which the energy levels are decorrelated. To do so, 
we make the approximation that the only effect of interactions is to compress the spectrum by 
an amount ((1 + g^/vfY — ((72/'^/)^))^^^, where g2 and 5^4 are typical energy interactions. This 
compression of the spectrum means that the effective mass increases if g2 is large compared 
to 5^4. This approximation corresponds to taking local interactions in the real space. Then, 
we see that the integrable modes under request are such as 




^/kl< 1 + - - - ksT. (28) 



These modes are integrable provided they lead to well-defined bosons. We say that a boson is 
well-defined provided a sufficient number of particule-hole excitations enter into the summation 
(|2^), that is if 2A — |g| > 27ra/L, where a is a dimensionalless coefficient, which counts the 
number of non-zero terms in (^^, and is thus a measure of integrability. We get that 

keT ~ (2A - 2na/L) ((1 + g./vf^ - {g2/vff) (29) 

The energy scale fc^T* thus increases as the cut-off A increases, and decreases if the interaction 
strength g2 increases. 

3 Level spacing statistics of a 2 D spinless fermion sys- 
tem in the A-restricted Hilbert space 

3.1 Bosonization of the bidimensional Fermi surface 



The bosonization of the Fermi liquid |jT9|, |20|, ^T[| involves a covering of the Fermi surface by 



spheres of radius A. Each small sphere is labelled by an integer a and one assumes that the 
Fermi surface is fiat at the scale A. Assuming that A ^ 27r/L, the commutation relations of 
the operators 

Pq,a = E ^"(k + q/2)ea(k - q/2)nq(k), (30) 
q 

with 

= C+_q/2Ck+q/2, (31) 

must include a Schwinger term, leading to the commutation relation 

[Pq,">Pj,a'] = '^a,a'WEC'"(k,q)K+H "^-h), (32) 

k 
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where the constraint Cn is 



C„(k, q) = e^ik + |)^,(k - |)0,(k' + |)^,(k' - |). (33) 



To simphfy the expression (|32D, one makes the assumption of a flat Fermi surface in each 
sphere of radius A. In order to specify this condition, one imposes that no bosonic excitations 
with an angle 6 > 7i/2 exist, where 6 denotes the angle between kf and q. The maximal angle 
theta is such as 

ta„(«-f)<|j, (34) 

which leads to 

A < (35) 

The equation (^) means that the curvature is negligeable in a sphere of radius A, and is 
compatible with the condition 

A > 2n/L (36) 

provided the linear length is large enough: L ^ 7c/kf. If A <^ 27r/L and |q| <^ A, the leading 
term in ( P^ is of the form [pq,a, Pq/ ^'1 = ^a,a' Ci{q.na) . We evaluate the corrections under 
the conditions A ^ 2tt/L and |q| -C A. Under these assumptions, the leading term represents 
the number of states in the parallelogram of figure |^. To obtain the corrections, one has to 
substract the number of states contained in the small shaded triangle. The number of states 
to be removed is approximately equal to 

^ (^) l^l^^os^si^^ = ^"(^•^°)-^^^^^~- (37) 
The commutation relations are thus of the form 

I q A rif^ I 

[Pq,a' Pq',a'] = ^a,a'^q,q' V^a(q-na) { 1 + 0{ 7— ^)}. (38) 



Following reference , we define 

aq(kf) = ^ - kf |)[nq(k)0(q.Vk) + n_q(k)0(-q.Vk)], (39) 

k 

and 

6q(kf) = (iVAV|q.Vk|)-i/X(kf), (40) 

where 0a is a smearing function such as (f)\ — > 5k,kf if A — > 0, and A^a is the local density of 
states: 

NA = ^T.\M\^-^f\?S{fi-e^). (41) 
^ k 
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Provided the curvature of the Fermi surface is neghgeable, that is provided condition (^) is 
satisfied, one can bosonize the Hamiltonian (0) to obtain 



^ = E E |q.Vk|6+(kf)6q(kf) + — 5: /k,k',qn_q(kK(k'). (42) 



kf q,q.Vk>0 k,k',q 



Since the interactions are quadratic in the boson operators, one can diagonahze them via a 
generahzed Bogohubov transformation. However, in numerical computations, one can only 
treat a Hilbert space of size 2000. This means that the cutoff A must be reduced drastically, 



as well as the Fermi wave vector kf. In particular, the conditions (35) and ( p6D are no longer 
valid and one cannot diagonalize the Hamiltonian as described above. Two scenarios are 
candidates for the appearance of chaos in the spectrum of the bidimensional Fermi liquid. The 
first scenario is the effect of curvature. In the limit L — +00, and for a cut-off independent 
of L, the condition ( P5| ) is no more valid, and one cannot solve the model by bosonization. 
We shall not study this effect in the present paper. The second scenario comes from the fact 
that, even with a flat Fermi surface, the condition (^) may be violated, so that the system 
is no more integrable by bosonization. We shall study the last type of effect in the rest of the 
paper. 

3.2 Level spacing statistics of a 2 D spinless fermion system in the 
A-restricted Hilbert space 

We now turn to the bidimensional case in the presence of a momentum cutoff. The Fermi 
sea in its fundamental state is pictured on the insert of figure |^. We treated a Fermi sea 
of 5 electrons for a total of 29 available quantum states. The Hamiltonian is given by the 
expression (0). The interaction term can be split into two terms as follows 

= Hl + Hl (43) 
^0 = -T7E/k,k+q,q^^k+q^^k (44) 

^ k,q 

^1 = TfE E /k,k',q4+qCi^'_qCk'Ck (45) 



^ q^Ok,kYk+q 



(46) 



The term is of the same nature as the diagonal Landau interaction between quasiparticules. 
The only difference is that 5nk represents occupation numbers of renormalized quasiparticules 
in the Landau theory, whereas is the number operator of bare fermions. The Hamiltonian 
made up of the kinetic term (^ plus the term (|4^ ) has already been studied in reference 



111 and leads to Poisson statistics in two dimensions. In the limit in which the bidimensional 
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bosonization procedure is applicable, it has been established in that the Hamiltonian 
H + Hi + H\ is integrable, namely that the effect of + H\ is to renormalize the free theory 
without breaking the integrability. We question whether this property is still valid for a system 
of electrons such as the one drawn on the insert of figure To see this, we diagonalize the 
Hamiltonian + + H\ for the electron system of the insert of figure ^, in a sector of 
fixed total momentum. As expected from bosonization and from R.P.A. theory p4[, we 



obtain a collective bound state which detaches itself from the particule-hole continuum. The 
energy of the collective mode is greater as the continuum energy for repulsive interactions, and 
lower for attractive interactions. Because of the small value of the momentum cutoff A, the 
energy width of the continuum is small compared to the energy of the bound state, so that the 
bound state renormalizes the level spacing statistics of the continuum towards small values, 
since the mean value of the level spacing distributions is rescaled to unity. If we suppress 
the bound state from the spectrum, we obtain level spacing statistics in good agreement with 
Poisson level spacing statistics, as shown on figure |^. It shows that whereas the cutoff A 
is very small and the number of electrons small too, the result we obtain is consistent with 
the result predicted in the limit in which the bosonization procedure is valid, with a large 
cutoff compared to 211 /L. However, there is no contradiction with what has been done in 
one dimension. If we assume that in two dimensions the statistics are poissonian provided 
a sufficient number of truncated bosons (|30D have a sufficient number of non zero terms in 
their linear combination, we obtain the same criterium as in one dimension, namely that the 
statistics are poissonian provided A ^ 27r/L. Using the cross-over scale we obtained from the 
diagonalization of interactions, we reach the conclusion that the level spacing statistics of 
the Hamiltonian + H]^ should be of the G.O.E. type. This is indeed the shown on 

figure [TO. We attribute the fact that the level spacing statistics of if° + + Hi are poissonian 
to the presence of only a small number of off-diagonal terms, which are not numerous enough 
to change significatively the statistics. 



4 Conclusion 

In this paper, we have studied the effect of the momentum cutoff on the level spacing statistics 
of interacting Fermi systems. We found that the presence of the cutoff could change drastically 
the level spacing statistics of a finite size system, namely to drive the level spacing statistics 
from a poissonian shape to G.O.E. statistics. Using temperature-dependant level spacing 
statistics and one dimensional models, we have shown that the system of electrons in the 
presence of a cutoff is integrable at low energy. As far as g2 interactions are concerned, we 
obtained G.O.E. level spacing statistics, which evolve to a Poisson level spacing statitics as 
the temperature decreases from +00. In the case, we could reach higher values of the cutoff 



13 



since the two branches of the dispersion relation are decoupled. The infinite temperature 
level spacing statistics were intermediate between Poisson and G.O.E., and were driven to 
a Poisson shape as the temperature decreases. The main feature of a system of spinless 
correlated fermions in the presence of a momentum cutoff in the /c-space is that, inspite of the 
loss of integrability there subsists nearly uncorrelated levels at low energy. This conclusion 
is analogous to the conclusions of reference ||25[ for another model. In two dimensions, one 
has to distinguish between two cases. The first case corresponds to off-diagonal interactions 
only. In this case, the level spacing statistics are G.O.E. statistics. It is clear that in this 
case, no switching on procedure can connect the excitations of the gas and the excitations 
of the interacting system, since the nature of the two spectra is different, which means that 
the system is not at the Landau fixed point. In the second case, diagonal interactions of 
the Landau type are taken into account. The Poisson level spacing statistics are restored 
because the off-diagonal matrix elements are not numerous enough, whereas their amplitude 
is comparable to the amplitude of the diagonal matrix elements. The Fermi liquid behaviour 
of the Hamiltonian in the presence of diagonal interactions is thus not destroyed. This result 
is to be compared with the fact that the Hamiltonian is diagonal in the limit in which the 
bidimensional bosonization procedure is valid. In spite of the small number of fermions and the 
small number of orbitals in our numerical computations, we obtain a result in good agreement 
with the bosonization of the Fermi surface theory. It should be stressed that the limit in 
which the bosonization theory is valid is not the thermodynamic limit in which A is fixed and 
L —>■ +00, because of the constraint (0). The problem of the thermodynamic limit, with the 
condition (^) violated remains open. 

The finite temperature level spacing statistics seem to be an appealing tool for the study 
of the spectrum of strongly correlated electronic systems. From a finite size study, one can 
characterize the integrability at low energy, and the cross-over temperature measures how far 
the low energy degrees of freedom are from being integrable. However, it is clear that the level 
spacing statistics retain only information about the the symetries, and give no information 
about the decay of the correlation functions, so that, in principle, we cannot solve entirely the 
problem whether the bidimensional t-J model is a Fermi liquid or not. Nonetheless, we can 
characterize the degree of integrability of the low energy degrees of freedom. A forthcoming 
paper shall be devoted to the case of finite size t-J models |]TD|. 

The author acknowledges B. Dougo for underlying the role of the momentum cutoff, and 
J.C Angles d'Auriac for help with algorithms. 
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Figure captions 



Figure 1: 

Modulus of the Fourier transform of tlie spectrum in tlie thermodynamic limit |/oo(t)| for 5 
bosons. Each pole corresponds to one or several definite boson excitations. We have chosen 
Vf = 1 and L = 2tt. A given boson g > generates a sequence of poles at r„^q = 21111/ q. The 
poles are denoted by their corresponding rational number n/q. 



Figure 2: 

Finite size effects of the modulus of the Fourier transform of the spectrum |/L(r)|. We have 
chosen Vf = 1 and L = 27c, A = 11. As shown on the figure, we can attribute boson quantum 
numbers to some peaks. We recognize the fractions 1/5, 1/4, 1/3,2/5,1/2,3/5,2/3,3/4, 4/5. 



Figure 3: 

Evolution of the energy levels in the presence of (72 interactions with a momentum cutoff. The 
parameters of the model are Vf = 1 and L = 2tt. The cutoff is chosen equal to 3, so that we 
have 3 fermions on each branch, for 6 available quantum states. The interactions are of the 
form g2q = 5'2oexp {—q/R). with R = 6. The momentum sector is 1, leading to 45 states in 
the Hilbert space. Because of the particule-hole symmetry, the spectrum is symmetric, but 
with no level crossings. 



Figure 4: 

Level spacing statistics of ID spinless fermions with g2 interactions and a cutoff. The param- 
eters are f j = 1, L = 27r. In order to eliminate the particule-hole symmetry present on figure 
^ we chose a non symmetric cutoff for particules and holes: the number of fermions on the 
right branch is 4, and the number of quantum states is 9. On the left branch, the number of 
fermions is 5 for 9 quantum states. The Hilbert space contains 1052 states in the sector of 
momentum 1. The level spacing statistics are found to be well fitted by G.O.E. statistics. 
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Figure 5: 

Level spacing statistics of ID spinless fermions with g2 interactions and a cutoff at finite 
temperature. The parameters are Vf = 1, L = 2tt. In order to ehminate the symmetry present 
on figure ^, we chose a non symmetric cutoff for particules and holes: the number of fermions 
on the right branch is 4, and the number of quantum states is 9. On the left branch, the 
number of fermions is 5 for 9 quantum states. The Hilbert space contains 1052 states in the 
sector of momentum 1. The level spacing statistics are plotted for inverse temperature equal 
to /? = 0, 0.1, 0.5 A cross-over is found from G.O.E. statistics as /5 = to Poisson statistics as 
(3 increases. 

Figure 6: 

Level spacing statistics of 1 D spinless fermions with interactions and with a momentum 
cutoff. The parameters are chosen such as = 1 and L = 2tt. The value of the momentum 
cutoff A is 12 and the total momentum is 24. The Hilbert space contains 1185 states. The 
statistics is found to be intermediate between a G.O.E. statistics and a Poisson statistics. 

Figure 7: 

Level spacing statistics of 1 D spinless fermions with interactions and with a momentum 
cutoff. The parameters are chosen such as = 1 and L = 2tt. The value of the momentum 
cutoff A is 12 and the total momentum is 24. The Hilbert space contains 1185 states. The 
statistics are plotted for the inverse temperature equal to /5 = and [3 = 0.3. The statistics 
for P = 0.3 is close to the Poisson statistics. 

Figure 8: 

Representation of a sphere at the Fermi surface. The commutation relation ( ^2]) is proportional 
to the number of states contained in the intersection of the parallelogram and the sphere n. 
The leading order term takes into account all the states in the parallelogram. 

Figure 9: 

Level spacing statistics of the Hamiltonian + Hq + Hi. The level spacing statistics are 
poissonian if one excludes the bound state, and renormalized towards s = in the presence 
of the bound state. The linear size of the box is 2tt. The kinetic term e{k) is quadratic, with 
a mass equal to 3. The interactions have the form fk,k',q = /°exp {—\q\/R), with R = 6. The 
insert represents the Fermi sea which was used. The Hilbert space contains 1042 states in the 
sector of momentum P = (27r/L)(2, 1). 

Figure 10: 

Level spacing statistics of the Hamiltonian + Hi with the bound state excluded. The 
number of zero level spacing is equal to 1114 , leaving only 186 non zero level spacings for the 
statistics, which explains the important fluctuations of the level spacing statistics which agree 
with the G.O.E. shape. The insert represents the Fermi sea which was used. The Hilbert 
space contains 1042 states in the sector of momentum P = (27r/L)(2, 1). 
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